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Abstract 


W'e  evaluate  the  mesoscale  vertical  velocity  induced  by  stress  changes  in  the  surface  layer  as 
a  function  of  the  size  of  the  rough  patch  in  relation  to  environmental  parameters.  The  nature 
of  the  flow  perturbation  strongly  depends  on  the  relation  between  the  width  of  the  rough  patch 
and  the  two  natural  scales  of  the  flow,  i.e.  the  inverse  inertia  wave  number  and  the  inverse  of  the 
Scorer  parameter.  When  the  width  of  the  rough  patch  is  comparable  to  the  inverse  inertia  wave 
number  or  larger,  the  atmospheric  perturbation  is  trapped,  the  vertical  scale  equals  the  depth 
of  the  stress  surface  layer,  and  the  horizontal  scale  equals  the  Rossby  radius.  When  the  width  of 
the  rough  patch  is  larger  than  the  inverse  of  the  Scorer  parameter,  but  smaller  than  the  inverse 
inertia  wave  number,  the  atmospheric  perturbation  is  a  hydrostatic  gravity  wave  with  a  vertical 
wave  number  equal  to  the  Scorer  parameter.  When  the  width  of  the  rough  patch  is  comparable 
to  the  inverse  of  the  Scorer  parameter,  the  atmospheric  perturbation  is  a  propagating  lee  wave 
with  a  vertical  wave  number  equal  to  the  Scorer  parameter.  When  the  ambient  flow  is  strong 
over  a  small  rough  patch,  the  flow  is  irrotational.  The  same  limitations,  inherent  to  the  linear 
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gravity  waves  excited  by  the  forcing  in  the  atmosphere  (e.g.  mountain  waves,  gravity  waves 
initiated  by  convection,  etc.),  apply  to  the  mesoscale  perturbation  induced  by  a  rough  patch 


1  Introduction 

When  an  air  mass  approaches  a  region  where  there  is  a  substantial  increase  of  the  shearing  stress  in 
the  surface  layer,  the  air  speed  decelerates  in  the  lower  layer.  The  resulting  horizontal  convergence  is 
associated  with  rising  motion.  Such  situations  are  typical,  for  example,  in  coastal  urban  areas  when 
onshore  flow  occurs.  However,  it  may  also  be  of  significance  in  inland  urban  areas,  because  of  the 
contrast  with  surrounding  agricultural  rural  areas  or  when  there  is  a  contrast  between  prairie  and 
wooded  areas.  The  developed  vertical  velocities  may  trigger,  under  supportive  synoptic  conditions, 
convective  clouds.  The  features  of  the  induced  vertical  velocity  may  also  be  of  importance  in 
dispersing  pollutant. 

In  this  paper  we  approach  the  problem  of  the  vertical  velocity  which  arises  because  of  horizontal 
inhoraogeneities  in  the  surface  shearing  stress  in  the  atmospheric  planetary  boundary  layer.  This 
study  is  an  extension  in  more  general  terms  of  a  previous  paper  (Dalu  et  al.  1988),  where  we  report 
on  the  waves  generated  by  a  change  in  surface  roughness  {CSR). 

Hunt  and  Simpson  (1982)  provide  an  excellent  review  of  the  studies  reported  by  that  time, 
concerning  flow  perturbations  induced  by  a  roughness  change.  Furthermore  Hunt  (Belcher  et  al, 
1990)  and  at  the  EUROMEC  1990  meeting  presented  solutions  concerning  the  flow  perturbation 
induced  by  roughness  changes  within  and  around  the  stress  layer.  Our  work  generally  agrees  with 
his  results,  however  it  is  less  detailed  in  the  structure  of  the  perturbation  within  the  stress  layer, 
because  we  are  more  concerned  with  the  flow  perturbation  in  the  free  atmosphere  above.  Additional 
studies  reporting  on  the  impact  of  a  sudden  change  in  the  surface  roughness  on  the  horizontal  flow' 
are  given  by  Pendergrass  and  Arya  (1984). 

Claussen  (1987)  computed,  using  a  model  simulation,  the  vertical  velocity  due  to  a  CSR.  How¬ 
ever,  the  computed  vertical  velocities  were,  in  general,  very  sensitive  to  the  horizontal  grid  resolu¬ 
tion.  which  must  often  be  reduced  to  several  hundred  meters  in  order  to  appropriately  resolve  the 
related  vertical  velocity.  Using  a  very  coarse  horizontal  resolution,  Vukovich  and  Dunn  (1978)  in 
their  numerical  model  simulation  of  the  St.  Louis  urban  area,  suggested  that  the  surface  rough- 
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ness  has  only  a  small  effect  on  the  circulation  for  the  wind  speeds  used  in  their  study.  Alestalo 
et  al.  (1985),  using  a  hydrostatic  two-dimensional  model  with  a  grid  interval  of  4  km  simulated 
the  airflow  in  the  Baltic  shore  region  of  Finland,  and  found  a  maximum  for  the  vertical  velocity 
of  order  of  1  cm  due  to  the  CSR.  They  attributed  the  reported  increase  of  precipitation  in 
that  area,  in  the  absence  of  thermal  forcing,  to  the  vertical  velocity  induced  by  the  CSR.  Pielke 
(1974)  evaluated  the  magnitude  of  vertical  velocity  caused  by  a  CSR  over  Florida  using  a  11  km 
horizontal  resolution  model.  Although  the  magnitudes  were  small  (as  0.1  cm  5“*),  it  was  concluded 
t;  at  shallow  warm-rain  clouds  over  the  southeast  coast  of  Florida  could  result  due  to  this  mecha¬ 
nism.  Finally,  Roeloffzen  et  al.  (1986)  presented  a  steady  state  model  calculation  of  secondary  flow- 
patterns  forced  by  a  CSR.  Adopting  a  neutral  boundary  layer  and  using  a  refined  grid  resolution, 
they  suggest  that  frictional  effects  involved  with  a  CSR  at  a  coast  line,  can  lead  to  a  secondary- 
circulation  on  the  mesoscale.  They  suggest  that  this  forcing  is  a  factor  in  the  observed  coastal 
frontogenesis  active  in  the  early  fall  along  the  coast  of  the  Netherlands. 


2  The  Governing  Equation  for  the  Linear  Problem 


If  we  assume  that  the  process  is  stationary,  two-dimensional  and  Boussinesq,  then  the  primitive 
equations  in  linear  form  can  be  reduced  to  a  Scorer  type  equation  for  the  vertical  velocity  in 
non-homogeneous  form: 


k^ 


+  (l2.k^)u.  =  .ilL 
pi 


(1) 


Because  of  the  linearization  the  perturbations  may  be  underestimated,  however,  since  the  solutions 
are  continuous,  there  are  no  limitations  due  to  grid-size.  For  a  derivation  of  equation  (1).  s^^ 
Eliassen  (1977).  The  hat  denotes  the  Fourier  transform  of  the  variable,  k  is  the  horizontal  wave 
number,  =  f/C  is  the  inertial  w-ave  number  (/  is  the  Coriolis  parameter,  U  is  the  ambient  flow- 
perpendicular  to  the  change  in  surface  roughness,  and  I',  is  its  shear),  r  is  the  resulting  shear 
stress,  and  /  is  the  Scorer  parameter: 
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where  N  is  the  Brunt-Vaisala  frequency  and  b  is  the  buoyancy  of  the  environment.  Equation  (1) 
can  rewritten  as: 


=  G\k)^ 

pU 


(3) 


with  i'^(k)  =  k^^ — ^  and  G^{k)  =  ^ 


F  -  kl 


k^  -  kl 


In  the  wave  number  region  where  u^(k)  <  0,  the  waves  are  trapped  around  the  perturbing  source 
within  an  e-folding  vertical  distance  equal  to  The  vertical  wave  number  po  for  the  trapped 
waves  is: 


k^o(k)  —  li  i^{k)\  =  1^1  yp — p  when  0  <  j^*)  <  fco  or  when  /  <  jfci  <  oc  (4) 

In  the  wave  number  region  where  I'^ik)  >  0,  the  waves  propagate  away  from  the  perturbing  source 
with  a  vertical  wave  number  equal  to  pi.  The  vertical  wave  number  pi  for  the  propagating  waves 
is: 


pi{k)  =  i/{k)  =  k 


P-k'i 
k^  -  kl 


when  ko  <  |A:1  <  /. 


(5) 


Note  The  theory  of  the  mesoscale  vertical  velocity  induced  by  a  rough  patch  is  formulated  within 
the  framework  of  the  well  established  gravity  wave  theory,  therefore  solutions  in  the  presence  of 
variable  Scorer  parameter,  shear  in  the  ambient  flow,  regions  of  neutral  stability,  and  critical  levels 
can  be  easily  treated,  because  the  related  mathematical  tools  are  already  in  the  literature. 


2.1  Green  Functions  and  Boundary  Condition 

The  advantage  of  writing  the  solution  in  terms  of  Green  functions  is  that  a  variety  of  different 
vertical  profiles  of  the  stress  could  be  easily  studied  through  a  simple  convolution  integral.  Using 
Green  function  theory  (Stakgold,  1979),  we  seek  for  the  solutions,  gik,z  -  r'),  associated  with  the 
governing  equation  (1)  for  a  point  source  forcing  S{x',z'),  which  satisfies  the  boundary  condition: 

g{k,z-z')  =  0  when  z  =  0  (6) 


Then  the  vertical  velocity  for  a  given  forcing  is: 

,,G^(k)T,Ak,z') 
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The  Green  function  for  the  upward  propagating  wave,  which  satisfies  the  radiation  condition  (Som- 
merfeld,  1912  and  1948)  and  the  boundary  condition,  is: 


gi{k,z-  z')  =  [expi(*:x  +  ^i(A:)  (z  ~  z'j)  -  exp  i(kx  +  fii(k)  1^  +  z'l 

The  second  term  is  the  mode  reflected  by  the  ground. 

Remark  For  verification,  we  derive  the  boundary  value  Green  function,  gg^{k,z): 


(8) 


d  . 


Ji'=0 


=  expi(kx  +  tii{k)  z) 


(9) 


which  is  the  Green  function  for  a  radiative  wave  in  the  mountain  wave  problem  (Smith,  1979). 
The  Green  function  for  the  upward  trapped  wave  is: 

goik,z  -  z')  =  [exp(jA:i  -  Mo(^)  M  -  x'l)  -  exp(iA-i  -  //o(t)  ^  +  2'])]  (10) 

The  second  term  is  the  mode  reflected  by  the  ground. 

Remark  For  verification,  we  derive  the  boundary  value  Green  function,  gg^{k.z): 


9Bcik.z)  =  -  I— 5up(A:,2  -  z') 


d  . 
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=  exp(jA:i  ~  po{k)  z) 


(11) 


which  is  the  Green  function  for  the  trapped  wave  in  the  mountain  problem  (Smith.  1979). 


3  Atmospheric  Response  to  Stress  Changes  in  the  Surface  Layer 

The  stress  has  the  same  direction  and  opposes  the  ambient  flow;  furthermore,  for  simplicity,  we 
assume  that  the  stress  decays  linearly  with  altitude  within  the  stress  layer: 


t{x,z)  =  Tq  Heih-  F(x) 


(12) 


where  tq  is  the  surface  shear  stress,  F{x)  is  its  horizontal  distribution,  and  He  is  the  Heaviside 
function.  The  mesoscale  perturbation  depends  mainly  on  the  intensity  of  the  stress  and  on  its 
depth,  and  weakly  on  its  vertical  distribution.  The  relation  between  the  surface  stress  tq,  the  wind 
intensity  U,  the  surface  drag  Co,  and  the  shear  velocity  u*  is  given  by  (Panofsky  and  Dutton. 
1984): 


Cd  =  |-  =  -^  =  0  ( - j 


(13) 
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Here  u*  is  the  shear  velocity,  k  is  the  von  Karman  constant  and  zo  is  the  surface  roughness.  The 
atmospheric  response,  to  a  horizontal  distribution  of  the  stress  (^(t)  is  assumed  to  be  an  even 
function),  is  given  by: 


u:(x.2)  =  /o,  +  /i  +  /( 
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The  variable  w  is  the  amplitude  of  the  perturbation  of  the  vertical  velocity.  Here 
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Ik 


-  ik2)(/2_/t2) 

The  Wfi^{k,z)  waves  are  trapped  around  the  top  of  the  stress  layer: 

^{exp(-/xo|2-ft|)-exp(-/io|2:  +  /j|)}  cos(fci) 

The  xv^^{k,  z)  waves  propagate  away  from  the  top  of  the  stress  layer: 

w^^  (k,  2)=  ^  {sin(/ii(z  -  A[  +  *1)  -  sin(/ij  \z  +  h\F  kx)} 

The  tilde  denotes  the  cosine  Fourier  transform: 

F{k)  =  I  dx  F{x)cos{kx)  =>  F{i)  =  —  j  dk  F(k)cos(kx) 

Jo  Jo 


(15) 


(16) 
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4  Vertical  Velocity  Excited  by  a  Bell  Shaped  Stress 


We  assume  that  the  drag  coefficient  is  Cp  =  3-10'^,  that  the  depth  of  the  stress  layer  is  /i  =  300  m. 
and  that  the  environment  parameters  have  the  following  values  N  =  10“*  s~*  and  /  =  lO”'*  s“’. 
We  evaluate  the  vertical  velocity  induced  by  a  rough  patch  with  a  horizontal  extension  a.  in  relation 
to  the  ambient  wave  numbers  ko  and  1.  A  bell  shape  distribution  of  the  stress  is  ideal  for  this  kind 
of  analysis  (as  shown  by  Queney  (1947)  and  by  Smith  (1979)  for  the  vertical  velocity  induced  by  a 
bell  shape  mountain): 


Tzzix.z)  = 


Tof'(z-h)  0^ 


/l2 


a*  + 


i  ^  ro6(z-h)ira  ,  ,, 

'}  =  - ^2 - ^exp(-a^); 


(18) 
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Here  6  is  the  Dirac  function.  From  eq.  (14)  the  vertical  velocity  is  given  by: 

fko 


u-(x,z)  =  /o,  +  /l  +  /( 


dk  Go(k)u.'^(k,z)  a  exp{~ak) 


j  dk  Gi(k)u'^^(k,z)  a  exp{-ak)  +  j  dl:  Go(fc)u'^(^,  2)  a  exp(-alr)|  (19) 


(a)  -  When  la  >>  koa  >>  1,  /qj  >>  I\  +  the  perturbation  is  horizontally  and 
vertically  trapped  inertia  wave. 

Due  to  the  exponential  decay  of  the  Fourier  transform  of  the  bell  function  for  increasing  values  of 
the  wave  number,  when  the  rough  patch  is  large  and  the  ambient  flow  is  very  weak,  the  contributions 
of  the  second  and  third  integrals  are  negligible  in  comparison  to  the  contribution  of  the  first  integral: 


w(x,  z)  Rs  /oj  =  J  ^  ~  ~  ^1^  cos(fcx) 


a  j 

[  (a -1- ^|2  + h|)2  - 

(a+^|2-h|)2-i2 

[  [(a  +  ^i2  +  /l|)2  +  x2j^ 

[(a  +  fi2-hl)2  +  x2]' 

(20) 


This  solution  represents  an  inertia  wave  which  is  horizontally  and  vertically  trapped,  as  shown  in 
Fig. la.  When  the  inertial  wave  number  is  large  (in  Fig.la  koa  =  3  and  a  =  100km)  the  air  particles 
are  displaced  upward  and  northward  within  a  Rossby  radius  distance  upstream  of  the  rough  patch, 
then  the  restoring  Coriolis  force  brings  them  back  to  the  previous  location  through  an  inertial  os¬ 
cillation.  The  maximum  intensity  of  the  perturbation  occurs  at  2  =  /i,  and  monotonically  decreases 
above  it  as  in  the  Ekman  solution. 


(b)  -  When  la  >>  koa  =  0(1),  /q,  >>  I\  +  I03]  the  perturbation  is  vertically  trapped 
inertia  wave. 

In  this  case,  the  second  integral  does  not  contribute  significantly  because  of  rapid  oscillations 
of  the  sine  argument;  the  third  integral  does  not  contribute  because  of  the  exponential  decay,  thus 


w{x,z) 


exp  -k 


/k  +  M  \ 


-  exp  -k 


cos(l:i)(21) 
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This  solution  represents  a  number  of  vertically  trapped  inertial  oscillations  (Fig.lb).  Again,  the 
perturbation  monotonical  decays  with  altitude  above  the  stress  layer. 


(c)  -  When  la  >>  1  >>  koa,  /j  >>  /o,  +  /oji  the  perturbation  is  a  hydrostatic  grav¬ 
ity  wave. 

In  this  case,  the  trapped  wave  contribution  is  negligible  in  comparison  to  the  contribution  of  a 
vertically  propagating  hydrostatic  wave; 

tL'(x.z)vl\  =  w-  I  dk  [sin(/|z  + Aj  +  Arr)  -  sin(/|c  -  h|  +  A:j)]exp(-cA:) 

2  Jo 

=  %TT- — jr  {[asin(/|’  +  h\)  +  zcos(/|r-|-  h|)]  -  [asin(/|c  -  hi)  -t-  zcos(/U  -  h])]}  (22) 

2  [a^  + 

The  hydrostatic  gravity  wave  is  shown  in  Fig.2a.  When  the  inverse  of  the  Scorer  parameter  is 
smaller  than  the  extension  of  the  rough  patch  (fa  =  3  and  a  =  1  km  in  Fig.2a),  the  perturbation 
has  a  wave  structure  with  a  vertical  wave  number  equal  to  /.  The  maximum  intensity  of  the  per¬ 
turbation  occurs  at  the  center  of  the  rough  patch  at  z  =  h. 


(d)  -  When  la  =  0(1)  >>  ftoo,  /j  -(- /qj  >  /o,;  the  perturbation  is  a  non-hydrosiatie 
lee  wave. 

For  this  situation,  the  contribution  of  the  first  integral  /o,  is  negligible,  the  vertical  velocity  is: 

u’(i,z)  s: /i -1- /oj  (23) 


However,  most  of  the  contribution  is  in  the  propagating  non-hydrostatic  wave; 

I  |sin( -  A:^|z  -f-  fij  -f-  kx)  -  sinfv^/^  -  A-^jz  -  kj  +  A'z)j 


'  “i/„” 


dk 


y/F^ 


1-2 


exp(-aA') 


TT  la 


u'— —  exp(-/a)  {[sin(/|z  +  /i|)  Jo{l{x  -  k  +  fi|))  +  cos(/|z  /i|)  Ho(l{x  -  jz  -I-  h\))] 
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-  [sin(/|z  -  h\)  Ml(x  -\z  -  h\))  +  cos{l\z  -  h\)  Ho{1(t  -  |z  -  h\))]} 

The  contribution  of  the  trapped  non-hydrostatic  wave  is: 

roo  /  [exp-(v/P^|z  +  h|)-exp-(v^P^|z- /i|)] 

/02  =  tl’- dk  — i - y/kiZrft - ^exp(-aA-)  cos(A:z) 
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« -u’^^exp(-/fl)  A'o(/x)  {exp(-/|2  +  Aj)  -  exp(-/|z  - /i|)}  when  |/i|  >>  1 


The  trapped  wave  decays  exponentially  with  the  distance  from  the  top  of  the  stress  layer,  therefore 
it  interferes  with  the  propagating  wave  only  at  ^  as  A,  while  weakening  the  propagating  wave  up¬ 
stream  and  strengthening  the  propagating  wave  downstream  [the  zero  order  Neumann  function  is 
even,  while  the  zero  order  Struve  function  is  odd,  with  the  same  asymptotic  behavior  (Abramowitz 
and  Stegun.  1972)],  thereby  producing  a  wake  of  secondary  cells  downstream  at  the  level  of  the  top 
of  the  stress  layer.  The  non-hydrostatic  lee  wave  is  shown  in  Fig. 2b,  for  la  =  I  and  a  =  1  km).  The 
maximum  intensity  of  the  perturbation  occurs  at  the  center  of  the  rough  patch  at  r  =  A. 

(e)  -  When  koa  <<  la  <<  I,  Iq^  >>  /q,  +  /i,  and  the  flow  is  irrotational 

When  ambient  flow  is  strong,  the  inertia  wave  number  ko  and  the  Scorer  parameter  /  are  small. 
Thus  the  contribution  of  the  first  and  the  second  integrals  are  negUgible,  only  the  third  integral 
contributes; 

a 

w{x,  z)  =  Iqj  =  w-  I  tfA  [exp(-(a -I- [z -I- A|)A)  -  exp(-(a -I- |r  -  A|)A)]cos(Ai) 

2  Jo 

=  _ (a  +  ]z-h])  1 

2  l[(a-f  |z  +  A|)2-|-ar2]  [(a  +  |z  -  A|)2  +  ^  ’ 

The  perturbation  is  horizontally  and  vertically  trapped,  as  shown  in  Fig. 2c  for  la  =  0.3  and  a  = 

1  km.  The  maximum  intensity  of  the  perturbation  occurs  at  the  center  of  the  rough  patch  at 

z  =  A. 

5  Conclusions 

We  have  shown  that  a  horizontal  change  in  surface  roughness  can  induce  substantial  vertical  velocity. 
The  vertical  velocity  can  be  in  the  form  of  propagating  waves  or  in  the  form  of  trapped  waves:  in 
both  cases  the  perturbation  can  be  physically  relevant,  since  the  maximum  is  placed  at  the  top  of 
the  stress  layer,  i.e.  in  the  region  where  it  is  important  to  have  positive  vertical  velocities  in  order 
to  trigger  cumulus  convection. 

The  nature  of  this  perturbations  depends  on  the  width  of  the  rough  patch  relative  to  natural 
scales  associated  with  the  magnitude  of  inertia  wave  number  and  Scorer  parameter.  The  vertical 
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scale  is  related  to  the  ambient  Scorer  parameter,  when  there  is  vertical  propagation.  The  horizontal 
scale  is  related  to  the  Rossby  radius  for  weak  ambient  flow  over  larger  rough  patches.  VV’hen  the 
rough  patch  is  small,  the  horizontal  scale  is  related  to  the  inverse  of  the  Scorer  parameter. 

The  theory  which  we  use  is  derived  from  mountain  wave  theory,  (Queney,  1947;  Eliassen,  1977; 
Smith,  1979). 
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List  of  Figures 


Figure  la  -  Contours  of  vertical  velocity  induced  by  a  bell  shaped  distributed  surface  stress,  when 
the  inverse  of  the  inertial  wave  number  is  smaller  than  the  width  of  the  rough  patch;  koa  =  3, 
a  =  100  km,  under  weak  flow  condition  {U  =3  m/s,  tl’  =  1  cm/s  and  Au’  =  0.1  cm/s).  Figure  lb  - 
When  the  ambient  flow  is  U  =  10  m/s,  then  koc  =  1,  ii’  =  10  cm/s.  and  Au’  =  1  cm/s. 

Dashed  lines  represent  negative  contour  lines. 

Figure  2a  -  Vertical  velocity  isolines  when  the  inverse  of  the  Scorer  parameter  is  smaller  than  the 
width  of  the  rough  patch,  la  =  Z,  a  =  I  km.  under  weak  flow  condition  (C  =  3  m/s.  u’  =  1  cm/s 
and  Au'  =  0.1  cm/s).  Figure  2b  -  When  the  ambient  flow  is  U  =  10  m/s,  then  la  =  1,  u-  =  10  cm/s 
and  All'  =  1  cm/s.  Figure  2c  -  When  the  ambient  flow  is  strong  (t/  =  30  m/s),  then  la  =  0.3. 
u’  =  100  cm/s  and  AiZ’  =  10  cm/s. 

Dashed  lines  represent  negative  contour  lines. 
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